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Abstract 

3G and 4G mobile are based on CDMA technology. In order to increase the effectiveness of 
CDMA receivers large amount of effort is invested to develop suitable multi-user detector tech- 
niques. However, at this moment there are only suboptimal solutions available because of the 
rather high complexity of optimal detectors. One of the possible receiver technologies can be the 
quantum assisted computing devices which allows high level parallelism in computation. The first 
commercial devices are estimated by 2004, which meets the advert of 3G and 4G systems. In this 
paper we introduce a novel quantum computation based Quantum Multi-user detection (QMUD) 
algorithm, employing simple Positive Operation Valued Measurement (POVM), which provides 
optimal solution. The proposed algorithm is robust to any kind of noise. 

Keywords: Multi-user detection, Positive Operation Valued Measurement, Quantum computing. 
Quantum Signal Processing 



1 Introduction 

The subscribers of next generation wireless systems will communicate simultaneously, sharing the 
same frequency band. All around the world in 3G mobile systems apply Direct Sequence - Code 
Division Multiple Access (DS-CDMA) promising due to its high capacity and inherent resistance to 
interference, hence it comes into the limelight in many communication systems. Another physical 
layer scheme. Orthogonal Frequency Division Access (OFDM), is also often used e.g. for Wireless 
LANs (WLAN) or HiperLAN, where the subscriber's signal is transmitted via a group of orthog- 
onal frequencies, providing Inter Channel Interference (ICI) exemption. Nevertheless due to the 
frequency selective property of the channel, in case of CDMA communication the orthogonality 
between user codes at the receiver is lost, which leads to performance degradation. Single-User 
detectors were overtaxed and showed rather poor performance even in multi-path environment [p. 
To overcome this problem, in recent years Multi-User Detection (MUD) has received considerable 
attention and become one of the most important signal processing task in wireless communication. 

Verdu has proven that the optimal solution is consistent with the optimization of a quadratic 
function, which yields in MLSE (Maximum-Likelihood Sequence Estimation) receiver. However, to 



find the optimum is a A'^F-hard problem as the number of users grows. Many authors proposed sub- 
optimal linear and nonlinear solutions such as Decorrelating Detector, MMSE (Minimum Mean 
Square Error) detector, Multistage Detector, Hoppfield neural network or Stochastic Hoppfield 
neural network |jl], |2|, |3|, Q| , and the references therein. One can find a comparison of the performance 
of the above mentioned algorithms in Q . 

Nonlinear sub-optimal solutions provide quite good performance, however, only asymptotically. 
Quantum computation based algorithms seem to be able to fill this long-felt gap. Beside the 
classical description, which we recently use, researchers in the early 20*'' century raised the idea 
of quantum theory, which nowadays becomes remarkable in coding theory, information theory and 
for signal processing 

Nowadays, every scientist applies classical computation, using sequential computers. Taking 
into account that Moore's law can not be held for the next ten years because silicon chip transistors 
reach atomic scale, therefore new technology is required. Intel, IBM and other companies invest 
large amount of research to develop devices based on quantum principle. Successful experiments 
share up that within 3-4 years quantum computation (QC) assisted devices will be available on 
the market as enabling technology for 3G and 4G systems. 

This paper is organized as follows: in Section |2| we shortly review the necessity of multi-user 
detection, whereas in Section |^ the applied quantum computation method is shown. In Section ^ 
we discuss the used system model. In Section ^ the novel MUD algorithm is introduced and finally 
we conclude our paper in Section |6[ 

2 Multi-User Detection 

One of the major attributes of CDMA systems is the multiple usage of a common frequency 
and time slot. Despite the interference caused by the multiple access property, the users can be 
distinguished by their codes. Let us investigate a DS-DCDMA system, where the z^^ symbol of 
the k^^ (A; = 1,2,... ,K) user is denoted by hk{z). Applying BPSK modulation, the output signal 
of the k*^ user, denoted by qk{t), is given as 

L 

Qkii) = \fE'k ^ bk{i)sk{t - zT), 

z=-L 

where Sk{t) and Ej. is the continuous signature signal and energy associated to the /c*'* user, T is 
the time period of one symbol, and (2L + 1) is the size of a block, respectively. For the sake of 
simplicity we assume one path propagation channel, so the channel distortion for the k^^ user is 
modeled by a simple attenuation factor hk{t) = a^- This model, however, can be easily applied to 
more sophisticated channel models, as well. 

The received signal is the sum of arriving signals plus a Gaussian noise component and thus 
can be written as follows: 

K 

r{t) =^hk{t) * Qkit) + n{t) = 

k=l 

K L 

= Y.Y1 VE'kakbkil)skit - IT) + n{t), (1) 

k=l l=-L 

where K is the number of users using the same band, n{t) is a white Gaussian noise with a constant 
Nq spectral density. In case of signatures limited to one symbol length and if the system is properly 
synchronized, then there is no intersymbol interference. Consequently, index z can then be omitted 
from ([l|), yielding: 

K 

r{t) = ^ y^kakbkSkit) + n{t),t G [0, T). 
k=l 



A conventional detector contains k = 1, . . . , K filters which are matched to the corresponding 
signature waveforms and channels and calculates the following decision variable: 



EkCik / r{t)sk{t)dt. 
Jo 



The traditional "single-user" detector (SUD) simply calculate the sign of expression (Q) yielding 
b'^^^ = sign{6fc}. This method results in poor performance, as bk contains not only the signal 
transmitted by the k^^ user but an interference term generated by the other users: 



bk = bkQkk + ; biQki + nk 




signal . ' noise 

multiple access interference 



where Qki is defined as follows: 

rT 



and Uk = Jq Sk{t)n{t)dt is a zero mean white Gaussian noise due to linear transformation. The 
output of the matched filter in vector form is: 

b = Rb + n, (4) 

where b = [61, 62, ... , bx]'^, b = [61, 62, • • • > bx]'^ and n = [ni, n2, . . . , nx]'^, whereas R = [gki, k = 
1,2, . . . , K I = 1,2, . . . , K] is a symmetric matrix with property of diagonal dominance {qzz > Qzj 
^j = l,2,...,K j^i). 

Based on the model introduced above, we derive now the optimal MUD. The MUD algorithm 
processes vector b, which is the output of the matched filter. To obtain optimal solution based 
on Bayesian decision rule one wants to chose the maximal probability binary sequence conditioned 
by the received data series. The optimal Bayesian detection reduces to the following minimization 
problem |]l|: 



Qki = \fE'k\fE'iakOii I Sk{t)si{t)dt, (3) 

JO 



h°P^ = min 

ye{-i,+i}^^' 



(b-Ry)^R-^(b-Ry) . (5) 



Unfortunately, the search for the global optimum of (|5|) usually proves to be rather tiresome, 
which prevents real time detection (its complexity by exhaustive search is 0{2^)). Therefore, our 
objective is to develop new, powerful detection technique, which paves the way toward real time 
MUD even in highly loaded system. 



3 Quantum Computation Theory 

Quantum theory is a mathematical model of a physical system. To describe such a model we need 
to specify the representation of a system. Every physical system can be characterized by means of 
its states in the Hilbert vector space over the complex numbers C. The vectors will be denoted as 
|(^)[|. The inner product (ijl^) maps the ordered pair of vectors to C with the properties |^ 

• Positivity: (^|V')> for IV") = 0, 

• Linearity: ((^|(a|^/;i)+b|V'2))=a(v3|V'i)+b(i/?|^2), 

• Skew symmetry: {(flip) = (iflip)* . 
^Say ket (p. 



In the classical information theory the smallest conveying information unit is the bit. The 
counterpart unit in quantum information is called the "quantum bit" the qubit. Its state can 
be described by means of the state \ip), tp = a|0) + where a,/3 G C refers to the complex 
probability amplitudes and |ap + 1/3|2 = 1 ^, g |. The expression \a\ denotes the probability that 
after measuring the qubit it can be found in |0) computational base, and shows the probability 
to be in computational base In more general description an A^-bit qregister \ip) is set up from 
qbits spanned by \i) i = . . . (M — 1) computational basis, where M = 2^ states can be stored in 
the same time |9|. 

A/-1 
1=0 

where M denotes the number of states and Vi / = 0, = 1, X] IV'jP — 1' respectively. 

It is worth mention, that a U transformation on a qregister is executed parallel on all M stored 
states, which is called quantum parallelization. To the irreversiblity of transformation, U must be 
unity = U'^ . The quantum registers can be set in a general state using quantum gates which 
can be represented by means of a unitary operation, described by a quadratic matrix. Applying 
four basic gates any states can be prepared Q. 

4 System Model 

In DS-CDMA systems an information bearing bit is encoded by means of a user specific code with 
length of the processing gain (PC)]!!. In case of uplink communication we assume perfect power 
control. In the receiver side it is not required synchronization between input signals and user 
specific codes, however, chip synchronization is necessary, that allows using multipath propagation 
channels. 

Since classical multi-user detection schemes only try to minimize the probability of error in 
noisy and high interference environment, they, even also optimal solutions, can commit an error. 
Actually, these classical approaches make compromise between computational complexity, proba- 
bility of error and time barrier required for efficient working. On the other hand, QMUD does not 
make error in detection, at most QMUD can not make a decision in certain cases, furthermore, it 
indicate us this symbol in order to correct in a higher layer. 

4.1 Representation of Possible Received Sequences in Qregisters 

We quantize every chip of the k^^ user's codeword in a qregister of length N^h, where the number 
representation is not significant at the evaluation of the received symbol, because with increasing 
the number of the alphabet, the precision of number notation grows exponentially. In our model 
we prepare for user k two quantum register \ip^) and |(^q) each corresponding to transmitted bit 
"1" and "0" with an overall length Nq = N^h ■ PG. It is important to notice that the effects of 
a multi-path channel and the additive noise are contained in the registers, moreover, the density 
function of the noise does not need to be known a-priori, however, the knowledge of noise can 
reduce the number of qubits in a qregister. This uncertainty may not influence the exact decision. 
Let V denotes a vector space spanned by \vi), i = 1 . . . orthonormal computational base states, 
where {vi\vj)=0 for \/i ^ j and {vi\vj)=l for \/i = j is hold. The number of stored states in quantum 
registers \ipi) or \ipq) is denoted with A^^i or Ngo, respectively. If the register \ip^) contains the 
desired state \vi), then 

(.fK> ^ { % (7) 

that fulfills the stipulation ^^^j^ =1- 



4.2 Preparation of quantum register states 



Due to the effect of multi-path propagation it is required to form any delayed version of chip 
sequences of user k. This operation can be made via the so called swap gate, which changes the 
position of two qubits in a register. In general, it can be seen as a quantum shift register. One 
can think, all the possible states should be computed before doing quantum multi-user detection. 
It is true, however, using classical sequential computers, this operation could take rather long 
time, whereas quantum computation exploits the quantum parallelism. Applying this feature a 
transformation on N states stored in a register can be done in one single step, that provides fast, 
efficient preparation of \ipi) and Iv^q)- 

5 Quantum Multi-User Detector 

The decision rule of classical multi-user detector becomes a measurement in quantum world. In 
our case we have to find out that the received and quantized signal vector of user k \r'^) = \vi) is 
either in the register Itp'f) or |(/9q) or both. In more mathematical description 

((^fc|r^) = i.e. f'l{i)=0. (8) 

Because of the multi-path propagation and the noise the same state \vi) could be found in both 
registers that makes the detection impossible. It shall be emphasized, however, that QMUD is able 
to recognize this event allowing higher layer protocols to perform error correction, hence it will 
never made false decision, as classical MUD algorithms (independently whether it is suboptimal 
or optimal) may do. On the other hand this can not be seen as feebleness of QMUD since the 
classical MUD is also unable to make proper decision in such a situation. The decision rules of 
QMUD are showed in Table ||. From now onward we only focus on \(Pi), the operations on |(^q) 
are analogous. 



Table 1: QMUD decision rule table 







decision 








no message was sent 







the bit "0" was sent 







the bit "1" was sent 


7^0 


7^0 


no decision is possible 



5.1 Evaluation of {ip^\r^) - The measurement 

The evaluation of {ip'Hr^) is not a trivial task as this is not an unitary operation, as discussed in 
Section ^. In the register \(pi) there is only one state \r^) = \vi) we are interested in. However, from 
measurement point of view the overall state of the quantum register being is state \ipi) can be re- 

garded as a qubit. This qubit can be written as a|0) + (995^ [r )|l), where a = yYlj=i jytil^iU)] = 
{ipi\vj). This qubit contains two states |?7i) = |0) and \r]2) = y^ ^^"^ |0) + ^ corresponding 
to the probability amplitude of \vi) is in the register or not. It can be simply proved that \r]i) and 
lr/2) are not unambiguously distinguishable, because {r]i\r]2) 7^ 0. 



However, we can extend the computational bases and apply the so called Positive Operation 
Valued Measurement (POVM-see Appendix). We introduce three positive operators 



El 



a|l)(l| 




a 



(9) 



Eo 



N. 



si 



-|0) 




(1| 



and 



(10) 



E.q 



I -El 



Eo 



:-/3(l- 




(11) 



where I is the identity matrix. The operator (11) provides 



(12) 



besides the first two POVM measurement operators in (^, ^) are orthogonal to \rii) and |r/2), 
respectively, making the probabilities of measuring Ei and E2 



P(Ei)|[^^^ = (7?i|EiIr?i) = 0, 

P(E2)|[^2^ = (7?2|E2|r?2)=0, 



(13) 



where -P(Ej) 1 1^ refers to the probability of the event the Ej was measured if |r/j) had been received. 
In other words, if our instrument indicates Ei, only the information corresponding to the state 7/2 
could be sent, otherwise if E2 is indicated the received state must be rji. It is appreciable that 
the scale of uncertainty arising from POVM measurement is a function of E3. It is important to 
emphasize that detecting E3 we do not any false detection. To reduce this effect the free variables 
a and (3 in E3 should be set to zero which makes it to identity matrix. Unfortunately, in that case, 
the resulting matrix becomes to a non-positive definite one. 

5.2 Setting the variables a and (5 

The operator E is positive if ((^lEjy?) > for any \ip). A positive definite matrix has the form 

E3 = ([vA, B){A, B\)^ where in our case A = \J'^ ~ ^ — y^l — a — /?^1 — > moreover, 
according to ([ll|) the product should satisfy 



AB 



13. 



IN. 



which leads to 



a 



1 



(14) 



(15) 



that makes E3 positive. 

We assume at the moment the symbols "1" and "0" are transmitted with equal probabilities, 
therefore it is worth choosing the measurement probabilities P(Ei)| j^^^ and P(E2)| |^^^ to be equal. 



Lemma 1. // the probabilities measurement P(Ei)||^^^ = P(E2)||^_^\| then a = (5 furthermore 
a = \. 

Proof. 

P(Ei)||^^^ = (7?2|Ei|r/2) = ^, 

and (16) 
= (m|E2|m) = ^ 

Substituting a = (3 m (|l5|) one gets a quadratic function with roots of ai = Ngi — \/ Nsi {Nsi — 1) 
and a2 = '^Ngi, where the latter one is impossible since probability can not become greater than 
1. Moreover, ai converges very fast to ^ as A'^si goes to infinity. □ 

However, as the length A^^i of a register grows the resulting probability of detection becomes 



always smaller, due to the small angle between |t/i) and |r/2). POVM [10| is typically used in such 
situations, where the measurement can not be repeated. In our case, however, in our system the 
content of the registers \(pi) and [(^q) are constant during detection, allowing multiple measurements 
or even parallelization of them, which makes P(Es) smaller at every step. 

Although, lemma 1. shows that a = /3 is an expedient choice, with additional thoughts -P(E3) 
can be further reduce. Therefore, we assert the next theorem. 

Theorem 1. Using appropriate different values for a and (3 one can double the probability of 
detection. 

Proof. We apply two measurements parallel, where 



and 



maxP(Ei)|i \ =^maxa 



ma^P(E2)||^^^. 

Focusing again on the former case, of course, P{Ei2)\r)i and also P{Ei^) become very small. Since 
the bounds of a probability variable x must satisfy < P{x) < 1, so the bounds of a, < a < A'^gi 
are known, as well. Issued from (15) a takes negative values if /? becomes greater than 1, and the 
same is held for the opposite case, respectively. The possible values of a and /3 are depicted in 
Figure ||, where two linear function of /3 can be seen according to the numerator and denominator 

of (HD. 

The maximum value for max a = 1 is also conceivable from Figure |l|, that makes P (E2)|j^_^^ = 

numerator (/3) 

denominator {/3) 

Figure 1: The boundaries of a 

and P'(Ei)||^^\j = which is 2 • P(Ei)\^^^y The same techniques can be applied for 

maxP(E2)||j^^\i , and it is enough for decision if one of the two measurements or both results 
El or E2. □ 



One can make a secure decision whether Ei or E2 or both is indicated as well as the effect of E3 
is reducible with repeated measurements. 

For a right quantum decision a Measurement Block (MB^) can be built up employing two 
POVM operations according to Theorem 1, as depicted in Figure ^, where k, i refers to the user 
and computation base (P), respectively. In general a detection can be made by a Decision logic 
operating on the following rule: 

1. inputi G {Ei,E2,E3}; 

2. if 3 inputi = Ei, i = 1, . . . ,k ^ out= Ei 



3. else if 3 inputi = E2, i = 1, 

4. else no decision. 



out= E2 



In our case the Decision logics and the Selector units in Figure § and ^ use the decision rule table 
to be found in Table according to decision rule described above. 



maxP(Ei^ 



1^2) 



POVM 
El, E2, E3 



inputi 



POVM 
El, E2, E3 



input 2 



maxP(E2) 



1^1) 



Decision logic 



I Measurement Block (MB) 

Figure 2: Measurement block 



It is remarkable that in case of Out = E3 the measurement block can be fed back i.e. the 
operation can be repeated, which is normally equivalent to switch such a block serial. 



Table 2: Measurement block decision rule table 



Input 1 


Input 2 


Decision 


Observation 


El 


El 


El 


the received symbol is surely in (/?i) 


El 


E3 


El 


the received symbol is surely in \ip\) 


E2 


E2 


E2 


the received symbol is surely not in 


E2 


E3 


E2 


the received symbol is surely not in 


E3 


E3 


E3 


no decision 











IVI 




Decision I. 
















IVI K*^ 






IVI 


















IVI Kf^ 




Decision II. 
















IVI Kf^ 






IVI k;^ 





Selection 



Figure 3: Receiver built from measurement blocks 

6 Conclusions 

In this paper we presented a quantum computation based multi-user detector algorithm, which 
involves the Positive Operation Valued Measurement. The new method utilizes one of the possible 
future receiver technologies of 3G and 4G mobile systems, the so called quantum assisted comput- 
ing. QMUD provides optimal detection in finite time and complexity when classical methods can 
achieve only suboptimal solutions. Our task is in the future to examine and underline the in this 
paper given theoretical results with some simulations. 
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Appendix 1: POVM- measurement 



POVM is a common used type of measurement, which provides a secure decision, however, does not 
care about the state after the measurement. The probabihty, notable as p{m) = {(p\ Mm^Mm \(p), 

Em 

where E^n is positive definite i.e. (<y9|Eni|<^) > and ^„ Em = 1 must be satisfied. One can 
construct a POVM with three elements/bases {Ei,E2,E3} in such a way, that in case of Ei or 
E2 unambiguous decision is possible between two occurrence. If E3 is indicated we can not decide, 
however, in worst case the error correction is handled in a higher layer protocol. 



